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P A ^ .1  -> 

In  this  study, Ve  have  applied  the  technique/ to  Nevada  Test  Site  (NTS) 
explosions.  The  mean  estimated  depth  Is  changed  from  approximately  50  km  to 
approximately  0 km  with  standard  deviations  of  30  km  for  a well-distributed 
5-statlon  network,  and  20  km  for  a 9-statlon  network. 
fi 

-pointy  out  that  the  technique  can  be  In  serious  error  If  deep  earth- 
quakes are  used  to  determine  residuals  for  shallow  explosions  In  a source 
area  where  the  earth  structure  between  the  earthquake  and  surface  Is  different 
from  that  Implied  by  the  travel-time  table  used. 


We  also  showQhat  there  Is  no  evidence  for  change  of  travel-time  residuals 
with  time  for  arrivals  from  NTS  at  RKON,  NPNT,  BUL,  ane  PRE.  There  Is,  however, 
evidence  that  significant  changes  In  residuals  are  correlated  with  location  at 
Pahute  Mesa  £ind  that  the  changes  may  be  due  to  Interactions  with  a deep  voj^anlc 
plug  under  Pahute  Mesa. 
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ABSTRACT 


Travel  time  residuals  may  be  obtained  from  a least-squares  location  program 
which  is  run  witn  the  depth  constrained  to  the  known  "true"  value.  Wlnen  these 
residuals  are  used  as  travel-time  corrections  in  the  same  program  run  depth- 
free,  nearby  events  are  located  with  smaller  errors  in  depth.  An  elaboration 
of  this  technique  has  been  denoted  the  SRST  (Source-Region-Station-Time) 
technique  by  K.  Veith. 

In  this  study  we  have  applied  the  technique  to  Nevada  Test  Site  (NTS) 
explosions.  The  mean  estimated  depth  is  changed  from  approximately  50  km  to 
approximately  0 km  with  standard  deviations  of  30  km  for  a well-distributed 
5-station  network,  and  20  km  for  a 9-station  network. 

We  point  out  that  the  technique  can  be  in  serious  error  if  deep  earth- 
quakes are  used  to  determine  residuals  for  shallow  explosions  in  a source 
area  where  the  earth  structure  between  the  earthquake  and  surface  is  different 
from  that  implied  by  the  travel-time  table  used. 

We  also  show  that  there  is  no  evidence  for  change  of  travel-tiint  residuals 
with  time  for  arrivals  from  NTS  at  RKON,  NPNT,  BUL,  ane  PRh.  ■"‘’o-p  is,  however, 
evidence  that  significant  changes  in  residuals  are  correlated  with  location  at 
Pahute  Mesa  and  that  the  changes  may  be  due  to  interactions  with  a deep  volcanic 
plug  under  Pahute  Mesa. 
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INTRODUCTION 


Positive  determination  of  a seismic  event  as  deep  suffices  to  classify 
it  as  an  earthquake,  while  if  it  can  be  definitely  established  that  the  event 
is  shallow  and  if  is  small,  then  one  may  classify  the  event  as  an 

explosion  with  a small  probability  of  error.  Thus  the  problem  of  depth  is  of 
crucial  importance  in  the  positive  discrimination  of  both  earthquakes  and 
explosions. 

The  most  generally  applicable  technique  for  depth  determination  is  by 
use  of  the  Geiger  location  technique  applied  to  the  P wave  arrival  times. 

This  method  is,  however,  inaccurate  for  shallow  depths  because  P waves 
propagated  to  teleseismic  distances  depart  the  source  at  almost  vertical 
incidence.  To  make  matters  worse  there  are  biases  in  depth  estimation 
resulting  from  the  use  of  incorrect  travel  time  tables.  See  for  example 
Flinn  (1965),  Evernden  (1969),  and  Chiburis  and  Ahner  (1970). 

One  approach  to  the  problem,  used  by  Evernden  (1969),  was  to  define  a 
better  regional  travel-time  distance  relation  by  use  of  explosions  of  known 
depth.  For  4 NTS  explosions  he  was  able  by  this  technique  to  reduce  the 
average  depth  from  39  to  12  kilometers.  Evernden  states  that  a similar 
result  was  achieved  using  earthquake  data  to  derive  an  accurate  travel-time 
table,  but  he  does  not  say  how  he  overcame  the  difficulties  of  unknown  depth 
and  origin  time  inherent  in  this  approach. 

In  a related  but  slightly  different  technique,  Evernden  used  improved 

travel-time  tables,  plus  P estimates  of  origin  time,  to  obtain  an  accurate 

8 

location  in  latitude,  longitude,  and  depth  for  the  Fallon  earthquake  which 
was  40  km  northeast  of  the  SHOAL  explosion.  Using  the  residuals  from  this 
location  for  Fallon  he  was  able  to  compute  a depth  of  -1.7  km  for  SHOAL  as 
contrasted  to  the  known  depth  of  +0.4  km. 

Another  technique  suggested  by  Evernden  was  to  use  the  residuals  from 
master  events  whose  depth— free  locations  gave  the  same  answer  as  pP.  However 
he  goes  on  to  say: 

"In  practice,  master  station  residuals  have  been  computed  against 
solutions  restrained  to  the  D(pP)  values  of  depth." 


Evernd.n  used  this  approach  1„  Ka»chatka-K„rlls.  and  found  good  agraa^ant 
atwaan  pP  and  Gaigar  depths  for  nao  sample  earthquakes.  He  also  found  that 
the  master  event  had  to  be  ulthln  l*-3"  of  the  event  of  Interest. 

Vaith  (1971.  1973.  1974.  1975)  made  a systematic  practice  of  what  the 

above  quotation  Indicates  Evarnden  felt  he  was  forced  to  do.  He  called  It  the 

Sourca-Ragion-Statlon-Tlma  (SRST)  technique.  Using  events  In  Kamchatka-Kurlls 

w th  good  pP  control.  Veith  constrained  the  events  to  the  pP  depth  and  then 

located  the  events  in  latltud,  and  longitude.  The  residuals  from  the  new 

location  ate  then  used  to  locate  different  nearby  events  depth-free.  The 

results  agree  well  with  pP  depths.  Veith  developed  a computational  technique 

Which  allows  the  t»thod  to  be  easily  applied  In  practice.  For  each  station 

the  residuals  for  a suite  of  events  with  good  pP  depths  covering  the 

Kamchatka-Kurll  region  were  contoured  by  fitting  a polynomial  to  the  measured 

values.  In  this  way  the  residual  tor  each  station  Is  allowed  to  vary  in  the 

proper  way  as  different  events  are  considered  or  as  the  event  location  shifts 

during  convergence  of  the  epicenter  estimation  process.  An  essential  aspect 

of  Velth's  approach  Is  also  to  use  distance  and  aslmuth-dependent  station 
corrections. 


Evernden  (1969),  however,  provided  an  indication  that  this  general 
approach  can  fail  in  critical  applications: 

"Jhere  are  conditions  under  which  the  use  of  earthquake  data 
con  rol  locations  of  explosions  does  not  work.  Thus,  station 

3l^£  d.a  earthquake  ^ea^ 

toitka  did  not  reduce  the  standard  deviation  of  Longshot  data 

e en  though  the  earthquake  was  only  about  1“  from  the  Longshot  ’ 
explosion  point.  In  other  words,  there  was  no  correlation  of 
Longshot  residuals  and  residuals  of  a nearby  shallow-focus 
Aleutian  earthquake.  The  epicenter  of  the  earthquake  was  south 
of  the  islands  in  the  general  area  of  the  trench.  Apparently 
the  marked  contrast  in  crustal  and  shallow  mantle  chaLcterisUcs 


We  may  note  however  that  these  same  proolems  must  have  existed  in 
Kamchatka-Kurils  where  plunging  oceanic  plates  distort  the  travel-times 
but  where  Evernden  and  Veith  had  little  trouble  in  obtaining  the 
proper  depth  for  earthquakes.  Could  not  the  problem  be  that  residuals  which 
are  proper  for  even  relatively  shallow  earthquakes  (e.g.  30-60  km),  may  be 


unsuitable  for  very  shallow  events?  Consider  the  case  of  a region  whose 
velocity-depth  characteristic  perfectly  matches  that  used  in  deriving  the 
travel-time  table,  except  that  the  surficial  AO  km  have  a low  velocity.  Then 
residuals  obtained  from  earthquakes  below  40  km  will  be  zero  and  will  there- 
fore be  incorrect  for  surface  events.  The  P wave  from  a surface  event  to  a 
distant  station  will  arrive  early  relative  to  the  arrival  at  a nearby  station, 
since  the  vertically  departing  P wave  spends  less  time  in  the  low  velocity 
surface  layer.  Thus  the  epicenter  of  a surface  event  will  appear  to  be  deep. 
As  we  shall  see,  the  upper  several  hundred  kilometers  of  the  earth  under  NTS 
appear  to  have  lower  velocities  than  those  implied  in  the  Herrin  68  travel- 
time tables,  so  that  NTS  events  are  located  too  deep.  Thus  if  deep  NTS 
earthquakes  existed  and  were  used  to  determine  SRST  corrections,  one  would 
have  the  embarrassing  situation  of  accurately  locating  all  earthquakes  but 
if  an  explosion  occurred,  one  would  calculate  a substantial  depth  for  it. 

Although  the  above  quote  from  Evernden  (1969)  might  seem  to  indicate 


that  he  was  aware  of  this  trap,  he  nonetheless  in  the  same  paper  presented 
a complete  picture  of  accurate  depth  estimation  in  Kamchatka  using  the 
technique  which  he  showed  to  fail  for  LONGSHOT.  He  apparently  concluded  that 
the  failure  was  due  to  lateral  changes  in  geology,  and  not  to  the  fact  that 
the  earthquake  and  explosion  were  at  different  depths. 

Since  Veith  found  that  near  Kamchatka  shallow  events  overlying  deep 
events  from  which  SRST  corrections  had  been  determined  were  located  too 
shallow,  we  presume  that  the  overlying  strata  had  greater  velocities  than 
those  below  relative  to  the  Herrin  68  Tables.  We  are  presently  studying 
a method  of  averaging  the  P and  pP  residuals  from  deep  earthquakes  to 
determine  residuals  which  would  be  appropriate  for  surface  events. 

The  main  purpose  of  this  report  is  to  confirm  Veith 's  procedures  with  a 
large  data  base  of  shallow  events,  i.e.  explosions  at  NTS.  In  the  course  of 
the  Investigation  it  became  clear  that  various  subsidiary  topics  could  be 
simultaneously  discusst i;  and  we  have  in  fact  investigated  the  variability 
of  travel-time  residuals  as  a function  of  time  and  space. 
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LOCATION  RESULTS 


As  a test  of  the  SRST  technique  discussed  in  the  Introduction  we  have 
used  tha  data  developed  by  Chiburis  and  Ahner  (1970).  We  have  constrained 
the  Pahute  Mesa  events  to  their  known  depths,  used  program  HYPO  to  perform 
the  Geiger  location  using  only  stations  with  A ^ 16“ , and  averaged  the 
residuals  for  each  station.  Stations  at  distances  greater  than  16°  were 
selected  in  order  to  model  the  problem  of  interest,  teleseismic  location. 
Table  1 gives  tne  mean  Pahute  Mesa  Herrin  6d  residual  for  each  station 
together  with  the  number  of  arrival  time  readings  and  the  standard  deviation 
of  the  population  of  residuals.  Figure  1 is  a map  of  the  depth  within  90° 
of  NTS.  The  station  numbers  from  Table  I have  been  plotted  next  to  one 
station  location.  Figure  2 is  a map  of  NTS  showing  the  general  size  of 
the  Pahute  Mesa  test  site.  All  of  the  events  in  this  report  have  very  few 
observations  in  the  southwest  teleseismic  quadrant.  Situations  such  as 
this  are  not  without  interest,  however,  since  they  frequently  arise  in  the 
practical  analysis  of  events  in  the  USSR  and  China. 

If  a signal  travels  through  a layer  of  thickness  d,  velocity  at 
incidence  angle  i^,  then  the  differential  travel  time  with  respect  to 
traveling  through  another  layer  of  velocity  is 


V.V  cos  i 
12  o 

Comparing  the  Basin  and  Range  structure  of  Masse  et  al.  (1972)  to  the  average 

world-wide  structure  implicit  in  the  Herrin  tables,  (Engdahl  et  al.,  1968)  we 

see  that  d=200  km,  Vj^=8.0  km/sec,  ^^=7.5  km/sec  are  reasonable  values.  Since 

any  small  region  such  as  NTS  may  have  a regional  correction,  we  may  assume  that 

the  variation  of  residuals  from  NTS  will  vary  as  1.7/cos  i^CA)  + constant.  The 

function  i (A)  is  taken  from  Richter  (1968),  Appendix  V,  who  assumed  a surface 
o 

velocity  of  6.34  kra/sec.  Application  of  Snells  law  shows  that  i^(A)  will 

change  insignificantly  with  respect  to  the  present  application  if  the  velocity 

structure  between  40  and  200  km  is  decreased  from  8.0  to  7.5  km/sec. 

The  function  1.7/cos  i (A)+constant  has  been  superimposed  on  the  residuals 

o 

from  Table  I plotted  in  Figure  3.  We  see  that  this  formula  accounts  for  the 
data  trends  rather  well.  The  overall  slope  is  in  general  agreement  with  the 
observations  and  the  substantial  cluster  of  large  positive  residuals  around  20 
degrees  seems  significant  since  it  is  at  just  this  point  that  the  theoretical 
curve  is  beginning  to  curve  sharply  up. 
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Mean  Paliute  Mesa  SRST's  for  Shots 
Constrained  to  Known  Depth  but 
with  Latitude  and  Longitude 
Unconstrained 


Standard  Deviation 
of  (ibse  rvat  ions 


Station  Number  of  Average 
Observations  Residual 


1 RK-ON 

2 KC-MU 

3 P(;-BC 
A FLO 

5 JE-IA 

6 OXF 

7 EU2AL 

8 SHA 

9 CPO 

10  AX2AL 

11  AAM 

12  WH2YK 

13  ATL 

14  BLA 

15  BE-FL 

16  SCP 

17  CMC 

18  CEO 

19  OGD 

20  COL 

21  LPS 

22  HN>ME 

23  SV3QB 

24  NP-NT 

25  CIE 

26  GDH 

27  SJC 

28  CAR 

29  NOR 

30  TRN 

31  KTG 

32  AKU 

33  NNA 

34  ARE 

35  KEV 

36  VAl. 

37  LPB 

38  ESK 

39  00-NW 

40  KON 

41  NUR 

42  COP 

43  PTO 

44  mat 

45  TOL 

46  STU 

47  GG-GR 

48  PEL 

49  HAL 
30  SHK 

51  SEO 

52  WES 

53  BEC 

54  GUA 

55  Sl-BC 

56  ANT 

57  TOR 

58  EM-MO 

59  KIP 


The  average  residuals  in  Table  I were  then  applied  to  the  eight  Pahute 
Mesa  events  which  were  used  to  derive  them;  and  to  fourteen  Yucca  Flat  events. 
These  events  are  listed  in  Table  II.  The  dates  and  origin  times  of  these 
events  and  of  several  others  to  be  discussed  in  the  next  section  are  given  in 
Table  IV.  Because  of  our  small  source  region,  station  corrections  are  not 
needed  as  a function  of  distance  and  azimuth  and  our  procedure  will  therefore 
give  results  equally  as  good,  in  this  case,  as  the  more  general  procedures  of 
Veith  (1971,  1973,  1975). 

From  Table  II  we  see  that  the  average  depth  of  Pahute  Mesa  events  was 
reduced  from  82  to  18  kilometers;  while  the  average  depth  of  Yucca  Flat 
events  was  reduced  from  64  to  6 kilometers.  The  LRSM  and  WWSSN  data  used 
is  poorly  distributed,  typically  for  A > 25°,  10  stations  NE,  10  SE,  1 SW, 

5 NW;  for  16°  < A < 2L°,  4 NE,  0 SE,  0 SW,  1 NW. 

All  these  events  had  a large  number  of  reporting  stations.  In  practice, 
however,  many  cases  of  interest  are  detected  only  at  5 or  6 stations.  We 
have  ised  the  COMMODORE,  GREELEY  and  TAN  arrival  time  data  to  simulate  this 
situation  for  some  cases  of  interest. 

We  required  first  a well-distributed  network  of  5 stations,,  consisting 
ot  one  in  the  distance  range  16°-25°  and  in  the  Northeast  quadrint;  and  four 
stations  beyond  A = 25°,  one  in  each  quadrant.  None  of  the  or, ter  four 
stations  could  be  closer  to  each  other  in  azimuth  than  30°,  Sets  of  stations 
meeting  these  criteria  vjere  selected  at  random  from  the  COMMODORE  data  using 
a Monte  Carlo  technique;  and  locations  were  obtained  without  the  use  of  SRST 
corrections  using  program  SHIFT  360. 

The  mean  resulting  depth  was  57  + 19  (standard  deviation  of  the  population) 
km.  Out  of  a total  of  50  trials  the  minimum  depth  was  16  km  and  the  maximum 
was  91.  When  SRST  corrections  were  applied  the  mean  depth  decreased  to  6 + 24 
km.  The  minimum  was  -46  km  and  the  maximum  was  +49  km.  When  the  random  SRST 
runs  were  per'ormed  with  readings  from  events  GREELEY  and  TAN  we  found  the 
results  given  in  Table  III. 

The  experiments  with  SRST's  for  all  three  explosions  were  repeated, 
allowing  4 more  stations  to  detect  beyond  25°.  Without  SRST  corrections 
the  mean  depth  for  COMMODORE  was  52  + 13.  The  minimum  was  27  km  and  the 


-13- 


TABLE  II 


Depths  Resulting  from  use  of  Mean  Pahute  Mesa  SRST's  on 
Pahute  Mesa  and  Yucca  Valley  Events  Together  with  the 
Estimated  Mean  and  Standard  Deviation  of  the  Sample  Population 


Number  of  Stations 
used  in  Location  Depth 

Event  Name  Calculation  Without  SRST  With  SRST 

Pahute  Events 


Rex 

17 

64 

6 

Duryea 

20 

80 

20 

Chartruse 

31 

92 

19 

Greeley 

51 

75 

13 

Scotch 

26 

67 

1 

Knickerbocker 

20 

100 

36 

Boxcar 

17 

68 

22 

Benham 

12 

107 

25 

82+16 

18+11 

Yucca 

Events 

Fore 

27 

36 

-17 

Klickitat 

20 

59 

34 

Turf 

21 

58 

0 

Wagtail 

30 

60 

13 

Cup 

20 

59 

16 

Buff 

28 

61 

2 

Dumont 

39 

84 

0 

Piledriver 

37 

59 

1 

Tan 

39 

70 

7 

Nash 

20 

65 

- 7 

Bourbon 

17 

80 

7 

Commodore 

38 

70 

3 

Auk 

14 

73 

34 

Corduroy 

38 

61 

- 3 

64+12 


6+14 


TABLE  III 


Mean  Depth  and  Standard  Deviations  of  the  Sample  Population; 
Minimum  and  Maximum  Depths  Out  of  50  Trial  Locations  using 
Observed  Data  in  Randomly  Selected,  Well-Distributed  Networks 
of  5 or  9 Stations;  with  and  without  SRST's 


Number 

Stations 

Base  Event 

Mean  Depth  and 
Standard  Deviation 
for  50  Trials,  km 

Minimum 
Depth,  km 

(50  Trials) 

Maximum 
Depth,  km 

(50  Trials) 

Without 

SRST 

5 

Commodore 

57  + 

19 

16 

91 

9 

Commodore 

52  + 

13 

27 

73 

With  SRST 

5 

Commodore 

6 + 

24 

-46 

+49 

9 

Commodore 

2 + 

16 

-22 

+27 

5 

Greeley 

15  + 

36 

-98 

+83 

9 

Greeley 

- 4 ± 

33 

-83 

+36 

5 

Tan 

-19  + 

41 

-87 

+44 

9 

Tan 

-16  + 

20 

-76 

+16 

5 

Average 

1 + 

33 

-77 

+58 

9 

Average 

- 6 + 

23 

-60 

+25 
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maximum  73  km.  With  SRST  corrections  for  the  average  of  all  three  explosions 
we  obtained  for  the  mean  and  standard  deviation  of  the  sample  population 
1 + 33  km.  for  5 stations  and  -6  + 23  km  for  9 stations. 

These  results  suggest  that,  even  with  SRST's,  depth  estimation  from 
teleseismic  P arrival  times  will  not  be  a reliable  discriminant  for  events 
with  estimated  depths  of  less  than,  say,  50  km  detected  at  less  than  10 
stations.  It  must  be  admitted,  however,  that  the  distribution  of  stations 
in  these  simulations  is  not  ideal,  and  that  a higher  percentage  of  stations 
with  A < 25®  might  help  substantially. 


t 

STABILITY  OF  RESIDUALS 

Table  IV  gives  arrival  times  for  NTS  events  at  several  stations,.  These 
data  can  be  used  to  investigate  the  question  of  stability  Oi.  travel~time 
residuals  as  a function  of  time.  The  possibility  that  residuals  might  vary 
as  a function  of  time  due  to  dilantancy  has  been  suggested  by  Wyss  (1973). 

We  decided  to  measure  the  time  of  arrival  of  the  first  maximum.  Illustrations 
of  signals  with  good  signal-to-noise  ratios  showing  the  picked  arrival  time 
are  given  in  Figure  4.  Since  all  the  signals  at  a single  station  were  in 
general  quite  smaller  (except  Pahute  and  Yucca  to  RKON  as  shown)  these  shapes 
were  used  as  a visual  "match  filter"  in  order  to  allow  the  analyst  to  more 
accurately  pick  the  arrival  time  for  weak  events.  We  have  plotted  times  for 
core  phases  which  because  of  their  small  angle  of  incidence  can  provide  good 
control  on  depth.  Both  PRE  recordings  in  Figure  4 show  precursors  which  we 
have  chosen  to  ignore  in  our  analysis,  picking  instead  the  sharp  arrival 
indicated. 

Figure  5 illustrates  the  travel-time  structure  of  Qamar  (1973)  in  the 
vicinity  of  BUL  and  PRE.  While  it  does  not  seem  to  provide  an  explanation 
for  the  precursors  mentioned  above,  which  arrive  too  long  before  the  main 
energy,  the  Figure  could  be  interpreted  to  explain  some  of  the  phasing 
after  the  arrivals  we  have  picked.  Figure  6 from  Sweetser  and  Blandford 
(1973)  shows  that  the  amplitudes  of  these  core  phases  are  half  a magnitude 
unit  greater  than  those  of  P phases  received  at  90°,  and  thus  are  etsily 
detected. 

In  Figure  7 we  have  plotted  the  Herrin-68  residual  for  events  in  the 
Yucca  Valley  to  RKON,  HPNT,  BUj.,  and  PRE.  Generally,  with  the  possible 

■ exception  of  late  1969,  there  appears  to  be  no  smooth  variation  with  time, 
and  there  seems  to  be  little  correlation  between  the  traces.  Similar  remarks 

■ seem  to  apply  to  the  plots  of  the  relative  residuals  in  Figure  8. 

Ine  situation  seems  substantially  different  however  when  we  plot  the 
Herrin-68  residuals  for  events  at  Pahute  Mesa.  Here  a very  strong  correlation 
is  apparent  between  stations  in  both  Figures  9 and  10.  A test  of  the 
significance  of  event  effects  in  the  relative  residuals  was  significant  at 
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Figure 


the  90%  level  although  not  at  95%.  Since  this  is  not  a station  effect,  nor 
an  effect  of  the  variation  of  regional  stress  with  time  (since  it  was  not 
I'bserved  for  the  Yucca  explosions)  we  conclude  that  it  may  have  something  to 
do  with  variation  of  the  geological  structure  under  Pahute  Mesa. 

In  Figures  11a, b,c  we  see  contour  map  interpretations  of  the  RKON  and 
BUL  residuals.  Although  it  is  not  clear  how  one  would  perform  a significance 
test  of  the  hypothesis  that  the  observations  are  correlated  with  location, 
it  seems  to  the  author  that  such  a correlation  exists.  The  existence  of  the 
pattern  for  relative  residuals  (Figure  11c)  shows  that  it  is  not  due  to 
erroneous  origin  times  or  depths  of  burial  since  both  of  these  effects  would 
be  expected  to  cancel  out  of  relative  residuals. 

It  is  worth  noting  that  for  every  station  the  variance  of  the  travel- 
time residuals  is  greater  for  Pahute  Mesa  than  for  Yucca  Flat  events.  This 
is  in  qualitative  agreement  with  the  arguments  of  Spence  (1973)  who  studied 
the  effect  of  an  hypothesized  Igneous  plug  under  Pahute  Mesa  by  calibration 
of  Pahute  Mesa  travel-time  residuals  against  the  apparently  simpler  patterns 
of  DUMONT,  a Yucca  Flat  event. 

A possible  explanation  for  the  fact  that  the  Pahute  residual  variation 
is  smaller  at  RKON  than  at  the  other  stations  is  that  rays  to  RKON  spend  a 
smaller  portion  of  their  time  in  the  high-velocity  volcanic  plug  since  they 
depart  at  a larger  angle  from  the  vertical  than  do  the  rays  to  more  distant 
stations.  Thus,  depending  on  how  close  the  event  was  to  the  center  of  the 
volcanic  plug,  the  residuals  to  RKON  and  other  stations  would  increase  or 
decrease  together;  but  with  a different  amplitude. 
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CONCLUSIONS  AND  POSSIBILITIES  FOR  FURTHER  RESEARCH 


Residuals  obtained  from  NTS  explosions  constrained  to  their  known  depth 
and  then  located  by  standard  Geiger  methods  may  be  used  successfully  to 
determine  the  depth  of  new  shallow  events.  This  represents  confirmation  of 
the  results  of  Veith  (1971,  1973,  1974,  1975)  obtained  with  earthquake  data. 

The  standard  deviation  of  20-30  km  for  5-9  station  location  of  individual 
events  suggests  that  depth  determined  using  this  technique  will  not  be  a 
reliable  discriminant  for  weak  events  for  depths  less  than  50  km. 

There  is  no  apparent  variation  of  travel-time  residuals  with  time  for 
compressional  waves  from  NTS  to  RKON,  NPNT,  BUL,  or  PRE.  All  of  these  stations 
are  in  essentially  aseismic  areas.  It  is  conceivable  that  residuals  to  a 
station  in  a seismic  area  such  as  MAT  in  Matsushiro,  Japan  would  e;  uibit  a 
variation  with  time  associated  with  the  dilatancy  effect,  Wyss  (1973). 

Although  one  cannot  reject  the  hypothesis  that  the  mean  Pahute  Mesa 
residuals  are  the  same  as  the  mean  Yucca  residuals  for  RKON,  NPNT,  BUL  or 
PRE,  there  is  a substantially  greater  variation  of  travel-time  residuals  at 
each  of  these  stations  for  events  at  Pahute  Mesa  than  for  events  at  Yucca. 

The  variation  appears  to  be  correlated  with  location  of  the  event,  and  is 
consistant  with  the  hypotnesis  of  Spence  (1973)  cf  more  complicated  structure 
under  Pahute  Mesa  than  under  Yucca  Flat.  The  standard  deviation  of  travel- 
time residuals,  even  over  as  small  a source  region  as  Pahute  Mesa  can  be  as 
large  as  0.4  seconds. 

With  respect  to  future  work,  we  are  presently  studying  the  possibility 
of  correcting  travel-time  residuals  computed  from  deep  earthquakes  to  the 
values  they  would  have  had  had  they  been  shallow  events.  The  approach  is  to 
average  P and  pP  residuals. 
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